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Exact stationary axially symmetric solutions of the 4-dimcnsional Einstein equations 
with co-rotating pressureless perfect fluid sources are studied. A particular solution 
with approximately flat rotation curve is discussed in some detail. We find that simple 
Newtonian arguments over-estimate the amount of matter needed to explain these curves 
by more than 30%. The crucial insight gained by this model is that the Newtonian 
approximation breaks down in an extended rotating region, even though it is valid 
locally everywhere. No conflict with solar system tests arises. 



1. Introduction 

The problem of the internal dynamics of galaxies is a very old one. It is of spe- 
cial interest because observed non-Keplerian (henceforth "flat" ) rotation curves ^ 
apparently contradict the simple Newtonian results. These observational facts are 
usually interpreted in favor of changing the matter model, i.e. introducing dark mat- 
ter that accounts for the difference between prediction and "experiment", cf. e.g.^1. 
However, one might also shift the burden to the gravitational theory itself. The 
most radical change in this line of arguments would result in avoiding dark matter 
altogether, as for example in MOdified Newtonian Dynamics^. So far these alter- 
natives typically are at least as unsatisfactory as the dark matter hypothesis itself, 
cf. e.g. ™. 

Taking a closer look at the usual approach, it is Newtonian Gravity rather than 
General Relativity (GR) that is applied to the situation of gravitational dynamics 
of galaxies. Because the velocity of the stars within a galaxy is small as compared 
to the speed of light and since the gravitational field is very weakfat least a couple 
of Schwarzschild radii away from a central galactic black hole ^ the Newtonian 
approximation seems to be sufficient. However, one has to be careful with such 
arguments. For instance, it is clear that corrections from GR are smaller in regions 
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where gravity is weaker - and yet, the GR distortion of the elliptic orbit of Mercury 
visually is most pronounced in the aphelion, and not in the perihelion, although 
gravity is weaker in the former. It is the main purpose of this paper to show that 
indeed the Newtonian approximation breaks down, even though it might be valid 
locally everywhere in the Galaxy. 

The idea to use GR rather than the Newtonian approximation to describe the 
internal dynamics of galaxies appears to go back to a relatively recent paper by 
Co ope rstock and Tieu^. Several errors were pointed out in subsequent publications, 
cf. and | Appendix B[ So it i s fair to say that Ref. ^ caused a lot of confusion, 
interest and critical comments I 7 | 8 | 9 | Actually, the approach of taking the non- 
linear character of GR more serious fits quite well to similar considerations on a 
cosmological scale 

Having a simple framework for the construction of exact solutions at our disposal 
helps to clarify these important issues. The prime goal of this paper is to provide 
and discuss such a framework and to point out where and why the Newtonian 
approximation breaks down. 

2. Derivation and solution of field equations 

2.1. Physical assumptions 

We do not discriminate between elliptic, spiral or irregular galaxies, but it is clear 
that some galaxies will obey the subsequent physical idealizations better than oth- 
ers. 

Order of magnitude estimates reveal that there are no relevant forces besides 
gravity in the description of intragalactic dynamics at the percent level. The matter 
content of a galaxy can be modeled very well by a pressureless perfect fluid. The 
approximation of vanishing pressure is applicable in the whole galaxy, aside from 
the very center where a galactic black hole may be located and a region close to 
the rotation axis where jets might be emitted. Additionally, stationarity and axial 
symmetry will be imposed, which is reasonable for many galaxies. This implies the 
existence of two Killing vectors and considerably simplifies the geometry. 

We are thus interested in axially symmetric stationary solutions of the Einstein 
equations (k = 87r, Gjv = c = 1) 

R ab ~\g ah R = KT ab , (1) 
with the energy momentum tensor 

T ab = pu a u b , (2) 
where p, u a are mass density and 4- velocity, respectively. 

2.2. Deriving the equations 

Theorem 7.1.1 of El (the premises of which apply to the case under consideration) 
guarantees that, without loss of generality, the line element can be brought into the 
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following form: 

ds 2 = -if (dt - iVd^) 2 + if "V d(j) 2 + n 2 {dr 2 + L dz 2 ) (3) 

All functions if, N, f2, L depend solely on the non-Killing coordinates r, z. The co- 
ordinate 4> is 27r-periodic. It should be noted that K = — f a £ Q is determined by the 
Killing norm of the timelike Killing vector £° = df. We are not interested in inside 
solutions of black holes here, whence £ a is timelike globally. 

We shall assume the perfect fluid velocity to be co-rotating, i.e., u a — (u°(r, z), 0, 0, 0). 
So actually u a is proportional to the Killing vector £ a = df. The conservation equa- 
tion V a T ab — decomposes into 

pu a V a u b = 0, V a (pu a )=Q. (4) 

If we suppose that p ^ then the first equation implies geodesicity of the 4- 
velocity while the second one expresses stationarity of the mass distribution. This 
seemingly mild assumption entails that the Killing vector £° has to be geodetic and 
hence if is constant. This is a crucial difference to the matter less case where no 
such restriction arises. Without loss of generality one may set if = 1 by a rescaling 
of the time coordinate in regions where p 7^ 0. Henceforth the vacuum case will 
always be considered as the limit of small p. Therefore, if = 1 is valid globally 
0. We emphasize that the presence of co-rotating pressureless perfect fluid sources 
simplifies the dynamics as compared to the vacuum case. 

In addition one may set L = 1. The proof as given e.g. in^^is valid only for 
the vacuum case, but one can show easily that the crucial condition R l t + Fft<j> = 
holds still for a pressureless perfect fluid (albeit it ceases to hold once pressure is 
switched on or a cosmological constant is added). Therefore, the line element ([3]) 
reduces to 

ds 2 =-(dt- Nd(f>) 2 + r 2 dcf> 2 + e u {dr 2 + dz 2 ) , (5) 

with the redefinition fl 2 = e v . From now on we shall exclusively refer to this adapted 
coordinate system. 

I nsertin g into the Einstein equations |TJ), J2]) establishes (cf. lAppendix A| and 
also™) 

rv z + N r N z = , (6) 
2rv r + N 2 - N 2 = , (7) 



1 

27" 2 

1 



Vrr + Vzz + 7h ( N r + N z) = > ( 8 ) 



N rr - -N r + N zz = , (9) 

r 

-j (A^ r 2 + N 2 ) = npe u . (10) 



a While this eliminates pure vacuum solutions of technical interest it retains the physically impor- 
tant solutions where "vacuum" is just an approximation to "very low energy density". 
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Remarkably, the only difference to the perturbative results in ^ is the factor e v on 
the right hand side in the last equation. 

Before solving these equations we would like to point out that the 3-velocity 
as seen by an asymptotic observer who is at rest with respect to the center of the 
galaxy is 

N(r,z) 

v ( r > z ) = — - — (11) 

provided r > N(r,z). This may be derived by a standard ADM split of the line- 
element 

ds 2 = -r 2 /{r 2 - N 2 )dt 2 + (r 2 - N 2 )[d(/) + N/{r 2 - N 2 )dt} 2 + e 1/ (dr 2 + dz 2 ) . 

(12) 

The standard decomposition of 9° into lapse and shift obtains 

d a t = 7 K + VE%) , (13) 

where n a = E? is the unit normal to the t — const, hypersurfaces and the shift 
vector is proportional to E^, as expected for angular rotation. The relativistic 7- 
factor is given by 7 = 1/ \/\ — V 2 with (fTT|) . So V is the velocity distribution of the 
(co-)rotating dust as seen by an asymptotic observer who is at rest with resp ect 
to the rotation axis. The results above can be found in the classic textbook D31 
(cf. also 01), but we have reviewed them in some detail here in order to make 
explicit the relevance of each of our assumptions, and also because the result (ITT]) 
for the velocity is crucial for our discussion. 



2.3. Solving the equations 

The strategy for solving this coupled system of partial differential equations (PDEs) 
is as follows: first, we find the general solution of the linear second order PDE © 
for N in terms of a spectral density. Then we use this solution to determine v from 
©-© up to an integration constant. Finally, we solve the last equation (fTO)) for 
p with N and v as input. To match with input from observations we must fix the 
spectral density (and thus the function N) accordingly. From (fTTj) we see that N 
is simply the velocity V times the coordinate r. Thus, the strategy is to take V as 
an input and to deduce the mass density p, rather than the other way around. 



2.3.1. Separation Ansatz and zero mode 

We implement now this strategy and focus in this Section on ©. The separation 
Ansatz N = lZ(r)Z(z) yields Z zz = kZ with the separation constant k £ E. Be- 
cause of reflection symmetry TV (r, z) — N(r, —z). Therefore, Z(z) = J dkA(k)Z(k, z) 
solely is composed of even modes Z(k, z) — cosh (y/kz). If k > then these modes 
diverge exponentially for \z\ —> 00. This is unphysical behavior. In^a positive value 
for k and fall off behavior at infinity was achieved by employing the non-smooth 
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modes Z{k, z) = e - ^^', which satisfy © for \z\ > only. Consequently, there are 
sources localized at z — 0. We show in | Appendix B| that either the weak energy 
condition is violated in the whole galactic plane or the pressure must be negative, so 
these sources are unphysical and contradict our assumption of vanishing pressure. 
Thus, we have no choice but to assume k = —A 2 with A <G Mq yielding the modes 
Z(X, z) — cos (Az). 

Let us address first the zero mode A = 0. Without loss of generality Z — 1 and 
1Z = Aq + Bqt 2 . Thus, the zero mode solution reads 

N = A a + B r 2 . (14) 

We consider first the simplest case, B = 0. Then the velocity (fTTj) is given by A /r 
and exceeds the speed of light for small r. While this is unpleasant, we note that 
the Einstein equations ([6|)- (fT0|) are invariant under the constant shift N — ► N — Aq. 
So the constant Ao may be absorbed into a redefinition of the time variable. Thus, 
in addition to the fall-off condition for \z\ — > oo exploited above and to the fall- 
off condition for r — + oo exploited below we have to impose a condition to fix 
Ao, e.g. by requiring N(0, 0) = 0. This will be discussed in more detail below; 
right now we simply fix Ao = 0. Then, if Bo ^ 0, the velocity profile V = Bor 
is linear, as expected from Newtonian gravity for a disk rotating with constant 
angular velocity. To get a first glimpse at possible differences between Newton and 
Einstein we consider the simplest case when only the zero mode is switched on, i.e., 
V oc r, and restrict ourselves to the galactic plane z = 0. In Newtonian gravity one 
obtains 

PN(r) = const. (15) 

However, general relativity behaves differently: plugging N oc r 2 into ||7J) yields 
v oc r 2 with negative proportionality constant. The unperturbed Einstein equation 
([T0|l establishes (a e R) 

p(r) oc e a2r2 . (16) 

Conceptually (although not physically) this is an important result: given as input 
the same velocity profile in the galactic plane, V oc r, Newtonian and general 
relativistic calculations yield different results for the mass density, even though 
locally the quantity p as given in (|16p can always be matched to the Newtonian 
in (|15[) by adjusting appropriately in v the integration constant a. Thus, only for a 
small region of spacetime the Newtonian picture is correct. However, the zero mode 
obviously is unphysical because for large values of r the velocity V would exceed 
the speed of light. Thus, from now on we set Bo = 0, i.e., the zero mode is dropped. 
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2.3.2. Generic modes 

In order to get the generic modes we use the rescaled variable x = Xr and obtain 
from the separation Ansatz above the ordinary differential equation 

K" - -K' - K = , (17) 

x 

where prime denotes differentiation with respect to x. It can be solved by standard 
methods. The result for N — HZ is given by 

oo 

N{r,z) = J dXcos(Xz)(rX)[A(X)K 1 (Xr)+B(X)h(Xr)] . (18) 
o 

The spectral densities A(X) and B(X) may be chosen according to the behavior 
required of N(r,z). The functions I\ and K\ are modified Bessel functions of the 
first and second kind, respectively (the latter is also known as Macdonald function). 
We note that I\ blows up exponentially for large values of r which is unphysical. 
Therefore, we set B{X) = 0. The function K\ falls off exponentially for large values 
of r and diverges like 1/r near r — 0. However, this divergence is compensated by 
a linear prefactor, so the integrand is well defined for any sufficiently regular A(X). 
We plot x ■ Ki(x) in Fig. |TJ The technical observation 




2 4 6 8 10 



Fig. 1. x-K 1 (x) 



J dxxKxix) cos (ax/r) = | +^3/2 ' ( 19 ) 


will be very useful below. 

To summarize, the general solution of © consistent with physical fall-off be- 
havior \z\,r — > 00 and with the absence of spurious sources in the galactic plane 
z = is given by 

00 

N(r, z) = A + J dX cos (Xz)(rX)A(X)K 1 (Xr) . (20) 
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The constant Aq has to be fixed by some additional boundary condition, e.g. by 
requiring a certain behavior of N as r — > 0. This can be done efficiently by setting 
first Aq — 0, designing a useful spectral density A(X) and finally shifting N by a 
suitable constant. 



2.4. Fixing the spectral density 

The choice of A(X) determines not only the radial velocity profile, but also the 
z-dependence of N. Thus, the behavior in z-dircction is determined uniquely by 
the behavior in r-direction and vice versa. This is good for the predictability of the 
model, because one may fit one kind of behavior and predict the other, in addition 
to predicting the matter density. Instead of fixing A(X) directly we found it useful 
to perform first a Fourier transformation, 

oo 

A(X) = | J dx C(x) cos {Xx) , (21) 
o 

where A(X) is determined in terms of a (Fourier) transformed spectral density C(x). 
This allows to represent V(r,0) as 

oo oo 

V(r 0)=r f dx ^ = - ( dx ° {ir) (22) 

v (r, uj r j ax + ^ )3/2 - ^ j ax + , ( zz) 



with some arbitrary function C'(x), which we may call again 'spectral density'. It is 
evident that for Coca; 2 the velocity profile is linear in r, while for C oc x the velocity 
profile is flat. Therefore, this purely technical step is very useful for matching with 
observational data. The function N(r > 0, z) reads 

oo 
2 /> 

N{r > 0, z) = r — i dxC{x)({{z + xf + r 2 y 3 / 2 + ((z - xf + r 2 )" 3 / 2 ) (23a) 
o 

For r = we obtain 

oo oo 

N(0, z > 0) = ^ J dX cos (Xz) J dxC(x) cos (Ax) = C{z) , (23b) 



provided the function C{z) is sufficiently well-behaved. The result for z < is 
obtained from reflection symmetry. 

For physical applications one still has to choose the spectral density C suitably. 
It is instructive to consider the case of constant C: the function N(r, z) turns out 
to be constant, which is problematic because for arbitrarily small values of r the 
function V(r, z) as defined in (|11[) grows without bound. In model building we 
have encountered this feature to be generic, rather than an artifact of a particular 
choice of C: if the velocity profile is taken from experimental data it turns out that 
for some values of z there will always be a region close to the axis r = where 
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V exceeds unity. This is a clear signal of the breakdown of our assumptions. For 
galaxies emitting jets this is not unexpected since jets emitted along the axis are not 
pressureless. Also the region around eventual central galactic black holes is located 
on the axis, so perhaps there is good (physical) reason for these regions not being 
described by a dust-like matter model. It is not the purpose of the current paper to 
further discuss these regions, where closed time like curves may arise, similar to the 
Godel universe (the appearance of such regions appears to be a generic feature of 
axially symmetric stationary solutions). Therefore, we cut out the domains where 
N(r, z) > r and blithely ignore them (cf. Fig. [2]). While this is clearly unsatisfactory, 
it appears to be not much worse than ignoring the problem with cusps in Newtonian 
cold dark matter models. 

3. A toy model for galaxies 

3.1. Input 

Experimental data imply ^ that for small values of r the velocity should increase 
linearly with r and for intermediate values of r it should yield an approximately 
flat profile. Asymptotically it should fall off like 1/r in order to match the Kerr 
solution H^l. Moreover, it is required that V(0, z) is regular for \z\ < ro, where the 
quantity ro is of the order of 1 kpc sa 3 • 10 54 , which allows to describe also the 
bulge region of the galaxy within our model. 

With these remarks in mind we propose a simple 3-parameter function for C (8 
is the step function): 

C(x) = (x - r ){6(x - r ) - 6{x - R)) + {R- r )9{x - R) (24) 

The three parameters, Vo,ro,R, have to be chosen appropriately. The quantity Vq 
is of the order of the velocity in the flat regime (Vq w 200 km/ s «7* 10~ 4 ). The 
bulge radius ro is set to 1 kpc and R to 100 kpc. Note that C is smooth besides 
the points x = ro,R. We stress that a smooth version of C may be chosen to 
improve the model (and also other refinements may be implemented according to 
phenomenological necessities), but we want to keep the model simple and exactly 
soluble for pedagogic reasons, while still being realistic enough to be applicable to 
galactic physics. 

3.2. Output 

From ([24]) and (|23ap we obtain 

N(r,z) = V (R - r ) + ^ £ (yj[z ± r ) 2 + r 2 - y/(z ± R) 2 + r 2 ) , (25) 

It may be checked that (f2"5")) indeed satisfies © up to some distributional contribu- 
tions at r = 0, \z\ > tq. The latter lie in the region we are not capable to describe. 
It is an open problem how to deal with that region along the axis (even in the 
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Newtonian approach), and we shall not address this issue here any further. How- 
ever, as we shall see now, these critical regions lie outside the toy galaxy, i.e., at 
values \z\ > ro, and are localized along the z-axis. Therefore, we are able to draw 
meaningful conclusions for galactic rotation curves. 




200 400 600 800 1000 



2 4 6 8 10 



Fig. 2. Contour plot of V (r, z) (N(r, z) > r in black region) 

Contour plots of the velocity are depicted in Fig. [2j The darker the region the 
higher the velocity. One unit corresponds to 1 kpc in these plots. The left figure 
shows that indeed the function V^(0, z) is regular until \z\ — 1 kpc, while the right 
figure reveals that on large scales the region we had to cut out stays close to the 
axis. The velocity profile in the galactic disk reads 

V(r, 0) = y (r - r + ^ r 2 + r 2 - V B 2 + r 2 ^j . (26) 

Various limits are V(r < r ,0) w rV /(2r Q ), V(r Q < r < i?,0) w V and 
V(R < r,0) w V R/r. Figure [3] depicts the ratio V(r,0)/V up to 20 kpc. The 




5 10 15 20 



Fig. 3. V(r, 0)/Vb plotted up to 20 fcpc 
velocity profile approximately behaves as required. Remarkably, ro determines also 
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the transition between the linear and the flat regime of the velocity profile, whereas 
the choice of R determines the transition between the flat and the asymptotic 
regime. Thus, within this model the thickness of the bulge may be predicted from 
the velocity profile. 

We have calculated the function v by solving ((6j)- dHJ and found that it is 
negligible up to the sub-permil level, besides the integration constant. This is 
mainly due to the fact that Vq /4 sa 10~ 7 <C 1. Thus, for most practical pur- 
poses v(r, z) ~ — ln/3 = const., where j3 is an integration constant. The density 
may now be derived from (|10| , 



Integrating the density over the white region in the right Fig. [2] yields about 10 11 
solar masses for the total mass, which is not unreasonable. The density within the 
galactic disk behaves as follows: for r — > is goes to a constant (PVq/tq), for 
r < r < R it falls off like 1/r 2 and for r > R it falls off like 1/r 6 . 

3.3. Comparison with Newton 

We shall now compare the mass density for a given velocity profile as predicted by 
simple Newtonian considerations with the GR mass density calculated above. We 
shall restrict ourselves to the galactic plane for simplicity. 

The Newtonian density for a given velocity profile may be derived naively from 
equating the kinetic energy per mass, V 2 /2, to the potential energy per mass, 
(47r J p N r 2 dr s )/r s . The quantity r s denotes the spherical radius r 2 = r 2 + z 2 , 
which coincides with the axial one in the galactic plane. Differentiation establishes 



The proportionality constant is known, but irrelevant as we are interested solely in 
the ratio p/ pn and p contains still an arbitrary multiplicative constant /3 which has 
to be fixed. 

In order to get a result which is robust with respect to slight changes in the 
spectral density C we assume an arbitrary function C with the following two prop- 
erties: 1. it leads to an approximately linear velocity profile for small r ("small" 
meaning, small as compared to 1 kpc, but not so small that one is at the scale of 
the galactic black hole, e.g. around 1 — 100 pc), 2. it leads to an approximately 
flat velocity profile for large r ("large" meaning, large as compared to sub-kpc but 
still smaller than the whole galaxy, e.g. around 1 — 20 kpc). For any given velocity 
profile V{r, 0) one may compare the Newtonian prediction with the one from GR 




(27) 



p N (r s ,0) = p N (r 7 0) oc 



V 2 + 2rVV 



(28) 



July 30, 2008 11:1 WSPC/ Guidelines Grumiller 



N on- Newtonian behavior in weak field general relativity for extended rotating sources 11 

B 



PN ' \ V 2 + 2rW 

Clearly, at each point in the galaxy we may achieve perfect agreement between 
Newtonian and GR predictions by choosing j3 appropriately. This reflects the lo- 
cal validity of the Newtonian approximation. However, it is important to realize 
that globally (3 can only be chosen once. Because usually one requires consistency 
between Newtonian and GR predictions in the region where luminous matter dom- 
inates, i.e., in the linear regime (V oc r), let us now fix [3 such that 



_P_ 

Pn 



linear 



I = •! J~. (30) 



Consequently, in the flat regime (V = const.) we obtain 



= P=\- (31) 

flat 4 



P_ 

Pn 

Thus, requiring consistency between Newton and GR in the linear regime neces- 
sarily leads to a considerable deviation in the flat regime: Newtonian estimates 
require 133% of the mass density as compared to GR calculations in order to get 
the same velocity profile. We reiterate the key observation: the integration constant 
j3 emerging in the GR solution can only be fixed once globally. 

We stress that these considerations are of relevance only for extended rotating 
sources. Thus, they play no role for solar system tests because there the main 
gravitational source is localized in the Sun. 



4. Conclusions 

Using a simple model for a galaxy we could show that the Newtonian approxima- 
tion over-estimates the amount of matter by about a third, whereas GR reduces 
the amount of dark matter needed to explain the flat rotation curves. This does 
not imply that exotic dark matter does not exist as insinuated by ®. After all, 
dark matter is a "500% effect" (meaning that there is about five times as much of 
dark matter in the universe as ordinary matter) and there are several independent 
measuremen ts a nd arguments which suggest the existence of dark matter, e.g. from 
microlensing ^3. But as astrophysical measurements reach unprecedented accuracy 
the unexpected GR correction in the 30% range predicted in our paper certainly 
will not be negligible, ft is emphasized that this prediction has been essentially 
independent from the particular choice of the spectral density C. 

It would be interesting to check the robustness of our results by relaxing some 
of our assumptions - in particular co-rotation - or (some of) the symmetry require- 
ments. 



b Obviously for the Keplerian profile V oc X/i/r the ratio l|29[ l becomes singular as a consequence 
of pjv = 0. 
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Appendix A. Ricci tensor 

In the adapted coordinate system ([5|) the components of the Ricci tensor read 

Rut, = |^ (-A^ r 2 + N 2 ) - r 2 {N rr - ~N r + N zz )^j (A.lb) 

Rtr = = Rtz (A.lc) 

^ ((r 2 + N 2 )(N 2 + N 2 Z ) + 2r 2 N(N rr - ±N r + N zz )) (A.ld) 

iV = = R+, (A.le) 

R rr = i (-u r r - v zz + ^N 2 + (A.lf) 

1/11 
2{r l 



R, 



R r* = ^i Z V * + — N rNz ) (A.lg) 



W 1 „ a 1 



Rzz = - \ ~Vrr - V ZZ + -jJV; - -V r j (A.lh) 

As usual, on the right hand side indices r and z denote partial derivatives with re- 
spect to the corresponding coordinate and repeated indices denote repeated partial 
derivatives, e.g. N rr := d r d r N(r, z). 



Appendix B. Localized exotic energy-momentum tensor 

In EH the expression (J% is a Bessel function of the first kind) 

N{r,z) = -^C(fc rl )\/fc„rJ!(\/fc„r)e-^H (B.l) 



July 30, 2008 11:1 WSPC/ Guidelines Grumiller 



N on- Newtonian behavior in weak field general relativity for extended rotating sources 13 

has been used, which is non-smooth at z — 0. According to the surface-layer formal- 
ism of Israel a localized energy-momentum contribution is generated by the jump 
in the extrinsic curvature across the surface z = 0, 

T l : h c = Kd(z)({K ab ] - h ab [K}) (B.2) 

where h ab denotes (the continuous) induced metric and [K ab ] — K^ b ~ K~ b the 
jump of the extrinsic curvature. Inserting into the definition (C n e denotes the Lie 
derivative along n c = E c z = e~ v / 2 d z ) 

K ab = ^£ n ch ab (B.3) 

with hob = -(dt-N d(j>) 2 ab + r 2 dct) 2 ab + e v dr 2 ab we find 

K ab = ^^(2(d< - Nd<f>)N z d0 + e v v z dr 2 ) ab . (B.4) 

The discontinuity [N z ] — 2j^ n C{k n )k n rJi{-\fk^r) inserted into ([6]) implies another 
discontinuity, [i/ z ] = — N r [N z ]/r. The jump of the extrinsic curvature 

[K ab ] = (2[iV 2 ](dt-7Vd0)d0 + e ,y K]dr 2 ) afc (B.5) 

yields the localized (effectively 1+1 dimensional) energy- momentum tensor 

^° b c = KS(z)^^N r [N z ]g ab := Pj~g ab , (B.6) 
with the 1+1 dimensional metric 

g ab = ~(dt -Nd^) 2 ab + r 2 dcf> 2 ab + ^-((dt -Ndtfi) dcj)) ab . (B.7) 

We note that both the 4- and 2-dimensional traces yield TrT^ c = 2Pj, whereas 
both the 4- and 2-dimensional contractions with the time- like Killing vector £ a yield 
the energy density 

Pi ■■= ffc?? = -Pi ■ (B.8) 

It saturates the strong energy condition, and additionally obeys the weak energy 
condition if N r [N z ] < at z — for all r. We may interpret the induced energy- 
momentum tensor as the one of a 1+1 dimensional perfect fluid with mass density 
pi and pressure Pj. So either the pressure is negative or the weak energy condition is 
violated. In any case we have exotic mattejfl. Moreover, one of our basic assumptions 
that led to the adapted line element J5]), namely the absence of pressure, clearly is 
violated by (|B~6]) . 



c One can decide between these two possibilities by taking a certain velocity profile in the galactic 
plane, inserting into calculating iV r and checking its sign. For monomials V = o?r a with 

a > —1 we get 7V r > 0. In addition, one needs an argument for the sign of [N z ]. 
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So the Cooperstock-Tieu model is unphysical and inconsistent, as pointed out in 
numerous comm ents . Where applicable, results of this appendix essential^ agree 
with the ones of 
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